
global data_dir "C:\Users\89384\Desktop"

//Figure 1 see Excel // 
 
//Figure 2 //
use "$data_dir\data.dta",clear
set more off

** Figure 2A
label define govern_c 0 "Non-Official" 1 "Official" 
label value govern_c govern_c
graph bar cig1 if retire==0 , over(govern_c) over (year) asyvars  ytitle(`"{stSerif:number/day}"',height(6)) title(`"{stSerif:A: Average number of cigarettes consumed}"') ylabel(0(2)10)

** Figure 2B
gen year1 = year if year == 2011 | year == 2013 | year == 2015
graph bar cig3 if retire==0, over(govern_c) over (year1) asyvars  ytitle(`"{stSerif:yuan/month}"',height(6)) title(`"{stSerif:B: Expenditure on cigarettes}"') ylabel(0(20)100)

** Figure 2C
graph bar f_liquor if retire==0, over(govern_c) over (year) asyvars  ytitle(`"{stSerif:number/month}"',height(6)) title(`"{stSerif:C: Average frequency of liquor drank}"') ylabel(0(5)15)

** Figure 2D
graph bar liquor if retire==0, over(govern_c) over (year) asyvars  ytitle(`"{stSerif:ml/month}"',height(6)) title(`"{stSerif:D: Quantity of liquor drank}"') ylabel(0(500)1500)

//Figure 3 //
use "$data_dir\data_PSM.dta",clear
set more off

gen test = year if antinew_governc == 1
bys ID:egen test_year = min(test)
gen distance = year - test_year

gen d_5 = 0
replace d_5 = 1 if distance == -5 | distance == -7
gen d_4 = 0
replace d_4 = 1 if distance == -4
gen d_3 = 0
replace d_3 = 1 if distance == -3
gen d0 = 0
replace d0 = 1 if distance == 0
gen d2 = 0
replace d2 = 2 if distance == 2
gen d3 = 0
replace d3 = 1 if distance == 3
gen d5 = 0
replace d5 = 1 if distance == 5

gen d_2 = 0

global cov age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area

** Figure 3A
reghdfe logcig1 governc d_5 d_4 d_3 d_2 d0 d2 d3 d5 $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
est sto r1
coefplot r1,keep (d_5 d_4 d_3 d_2 d0 d2 d3 d5) order(d_5 d_4 d_3 d_2 d0 d2 d3 d5) vertical recast(connect)  yline(0) omit scheme(plotplain) lc(gs6) levels(90) ytitle("LnNcig") ciopts(recast(rcap)) rename(d_5 = ≤-4 d_4=-3 d_3=-2 d_2=-1 d0=0 d2=1 d3=2 d5=3)  xline(4) note("Notes: 90%CI.")

** Figure 3B
reghdfe logcig3 governc d_5 d_4 d_3 d_2 d0 d2 d3 d5 $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
est sto r1
coefplot r1,keep (d_5 d_4 d_3 d_2 d0 d2 d3 d5) order(d_5 d_4 d_3 d_2 d0 d2 d3 d5) vertical recast(connect)  yline(0) omit scheme(plotplain) lc(gs6) levels(90) ytitle("LnEcig") ciopts(recast(rcap)) rename(d_5 = ≤-4 d_4=-3 d_3=-2 d_2=-1 d0=0 d2=1 d3=2 d5=3)  xline(4) note("Notes: 90%CI.")

** Figure 3C
reghdfe f_liquor governc d_5 d_4 d_3 d_2 d0 d2 d3 d5 $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
est sto r1
coefplot r1,keep (d_5 d_4 d_3 d_2 d0 d2 d3 d5) order(d_5 d_4 d_3 d_2 d0 d2 d3 d5) vertical recast(connect)  yline(0) omit scheme(plotplain) lc(gs6) levels(90) ytitle("Fliquor") ciopts(recast(rcap)) rename(d_5 = ≤-4 d_4=-3 d_3=-2 d_2=-1 d0=0 d2=1 d3=2 d5=3)  xline(4) note("Notes: 90%CI.")

** Figure 3D
reghdfe logliquor governc d_5 d_4 d_3 d_2 d0 d2 d3 d5 $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
est sto r1
coefplot r1,keep (d_5 d_4 d_3 d_2 d0 d2 d3 d5) order(d_5 d_4 d_3 d_2 d0 d2 d3 d5) vertical recast(connect)  yline(0) omit scheme(plotplain) lc(gs6) levels(90) ytitle("LnQliquor") ciopts(recast(rcap)) rename(d_5 = ≤-4 d_4=-3 d_3=-2 d_2=-1 d0=0 d2=1 d3=2 d5=3)  xline(4) note("Notes: 90%CI.")

//Table 1 //
use "$data_dir\data.dta",clear
set more off

** For the All sample
su cig1 cig3 f_liquor liquor anti official official_s governc anti_new age edu1 marrital gender health_condition wage_transfer fami_size area if retire == 0

** For the Difference between treated and untreated sample
ttable2 cig1 cig3 f_liquor liquor anti official official_s governc anti_new age edu1 marrital gender health_condition wage_transfer fami_size area if retire == 0,by(governc)

ttest cig1 if retire == 0,by(governc)
ttest cig3 if retire == 0,by(governc)
ttest f_liquor if retire == 0,by(governc)
ttest liquor if retire == 0,by(governc)
ttest anti if retire == 0,by(governc)
ttest official if retire == 0,by(governc)
ttest official_s if retire == 0,by(governc)
ttest governc if retire == 0,by(governc)
ttest anti_new if retire == 0,by(governc)
ttest age if retire == 0,by(governc)
ttest edu1 if retire == 0,by(governc)
ttest marrital if retire == 0,by(governc)
ttest gender if retire == 0,by(governc)
ttest health_condition if retire == 0,by(governc)
ttest wage_transfer if retire == 0,by(governc)
ttest fami_size if retire == 0,by(governc)
ttest area if retire == 0,by(governc)

//Table 2 //
use "$data_dir\data.dta",clear
set more off
global cov age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area 

   reghdfe logcig1 governc anti_governc $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m1
   reghdfe logcig3 governc anti_governc $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m2
   reghdfe f_liquor governc anti_governc $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m3
   reghdfe logliquor governc anti_governc $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m4

    esttab m1 m2 m3 m4 using "Table_2.csv",append order(governc anti_governc) legend starlevels( * 0.10 ** 0.05 *** 0.010) stats(N r2, star(r2) fmt(0 3) labels("Observations" "R-sq")) b(%9.3f) se(%9.3f) nogaps keep(governc anti_governc age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area _cons)


//Table 3 //
use "$data_dir\data_PSM.dta",clear
set more off

global cov age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area 

	reghdfe logcig1 governc anti_governc $cov if retire==0 [pw=w_ate], absorb(year communityID year#city_new) vce(cluster communityID)
	est store m1
	reghdfe logcig3 governc anti_governc $cov if retire==0 [pw=w_ate], absorb(year communityID year#city_new) vce(cluster communityID)
	est store m2
	reghdfe f_liquor governc anti_governc $cov if retire==0 [pw=w_ate], absorb(year communityID year#city_new) vce(cluster communityID)
	est store m3
	reghdfe logliquor governc anti_governc $cov if retire==0 [pw=w_ate], absorb(year communityID year#city_new) vce(cluster communityID)
	est store m4

	reghdfe logcig1 governc antinew_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
	est store m5
	reghdfe logcig3 governc antinew_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
	est store m6
	reghdfe f_liquor governc antinew_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
	est store m7
	reghdfe logliquor governc antinew_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
	est store m8
	
	esttab m1 m2 m3 m4 m5 m6 m7 m8 using "Table_3.csv",append order(governc anti_governc antinew_governc) legend starlevels( * 0.10 ** 0.05 *** 0.010) stats(N r2, star(r2) fmt(0 3) labels("Observations" "R-sq")) b(%9.3f) se(%9.3f) nogaps keep(governc anti_governc antinew_governc age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area _cons)


//Table 4 //
** DID estimations
use "$data_dir\data.dta",clear
set more off

global cov age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area 

	reghdfe logcig1 governc official_governc $cov if retire==0 , absorb(year communityID year#city_new) vce(cluster communityID)
	est store m1
	reghdfe logcig3 governc official_governc $cov if retire==0 , absorb(year communityID year#city_new) vce(cluster communityID)
	est store m2
	reghdfe f_liquor governc official_governc $cov if retire==0 , absorb(year communityID year#city_new) vce(cluster communityID)
	est store m3
	reghdfe logliquor governc official_governc $cov if retire==0 , absorb(year communityID year#city_new) vce(cluster communityID)
	est store m4

	reghdfe logcig1 governc official_s_governc $cov if retire==0 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m5
	reghdfe logcig3 governc official_s_governc $cov if retire==0 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m6
	reghdfe f_liquor governc official_s_governc $cov if retire==0 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m7
	reghdfe logliquor governc official_s_governc $cov if retire==0 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m8
	
	esttab m1 m2 m3 m4 m5 m6 m7 m8 using "Table_4_1.csv",append order(governc official_governc official_s_governc) legend starlevels( * 0.10 ** 0.05 *** 0.010) stats(N r2, star(r2) fmt(0 3) labels("Observations" "R-sq")) b(%9.3f) se(%9.3f) nogaps keep(governc official_governc official_s_governc age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area _cons)

** PSM DID estimations
use "$data_dir\data_PSM.dta",clear
set more off

global cov age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area 

	reghdfe logcig1 governc official_governc $cov if retire==0 [pw=w_ate], absorb(year communityID year#city_new) vce(cluster communityID)
	est store m1
	reghdfe logcig3 governc official_governc $cov if retire==0 [pw=w_ate], absorb(year communityID year#city_new) vce(cluster communityID)
	est store m2
	reghdfe f_liquor governc official_governc $cov if retire==0 [pw=w_ate], absorb(year communityID year#city_new) vce(cluster communityID)
	est store m3
	reghdfe logliquor governc official_governc $cov if retire==0 [pw=w_ate], absorb(year communityID year#city_new) vce(cluster communityID)
	est store m4

	reghdfe logcig1 governc official_s_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
	est store m5
	reghdfe logcig3 governc official_s_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
	est store m6
	reghdfe f_liquor governc official_s_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
	est store m7
	reghdfe logliquor governc official_s_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
	est store m8
	
	esttab m1 m2 m3 m4 m5 m6 m7 m8 using "Table_4_2.csv",append order(governc official_governc official_s_governc) legend starlevels( * 0.10 ** 0.05 *** 0.010) stats(N r2, star(r2) fmt(0 3) labels("Observations" "R-sq")) b(%9.3f) se(%9.3f) nogaps keep(governc official_governc official_s_governc age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area _cons)

//Table 5 //
use "$data_dir\data_PSM.dta",clear
set more off	
	
set seed 12345
set sortseed 12345
duplicates drop householdID year,force

reghdfe log_afood_overall governc anti_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
est store m1
reghdfe log_afood_overall governc official_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
est store m2
reghdfe log_afood_overall governc official_s_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
est store m3
reghdfe log_afood_overall governc antinew_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
est store m4

esttab m1 m2 m3 m4 using "Table_5.csv",append order(governc anti_governc official_governc official_s_governc antinew_governc) legend starlevels( * 0.10 ** 0.05 *** 0.010) stats(N r2, star(r2) fmt(0 3) labels("Observations" "R-sq")) b(%9.3f) se(%9.3f) nogaps keep(governc antinew_governc anti_governc official_governc official_s_governc age age2 edu1 marrital gender log_wage_transfer fami_size area _cons)
	
//Table 6 //
use "$data_dir\data_PSM.dta",clear
set more off		
	
gen need_control = 1 if smoke_control == 4
replace need_control = 0 if need_control == .

gen anti_smoke = anti*need_control
gen governc_smoke = governc*need_control
gen anti_governc_smoke = anti*governc*need_control

reghdfe logcig1 governc need_control anti_governc anti_governc_smoke $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
est store m1
reghdfe logcig3 governc need_control anti_governc anti_governc_smoke $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
est store m2
reghdfe f_liquor governc need_control anti_governc anti_governc_smoke $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
est store m3
reghdfe logliquor governc need_control anti_governc anti_governc_smoke $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)
est store m4
esttab m1 m2 m3 m4 using "Table_6.csv",append order(governc need_control anti_governc anti_governc_smoke) legend starlevels( * 0.10 ** 0.05 *** 0.010) stats(N r2, star(r2) fmt(0 3) labels("Observations" "R-sq")) b(%9.3f) se(%9.3f) nogaps keep(governc need_control anti_governc anti_governc_smoke age age2 edu1 marrital gender log_wage_transfer fami_size area _cons)	
	
	
//Table 7 //
use "$data_dir\data.dta",clear
set more off

global cov age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area

** Number of cifarettes
areg logcig1 governc anti_governc $cov i.year i.year#city_new if retire==0 , absorb(communityID) vce(cluster communityID)
psacalc delta anti_governc, rmax(0.409)
psacalc delta anti_governc, rmax(0.5)
psacalc delta anti_governc, rmax(0.6)	

** Expenditure of cigarettes
areg logcig3 governc anti_governc $cov i.year i.year#city_new if retire==0 , absorb(communityID) vce(cluster communityID)
psacalc delta anti_governc, rmax(0.434)
psacalc delta anti_governc, rmax(0.5)
psacalc delta anti_governc, rmax(0.6)

** Frequency of liquor
areg f_liquor governc anti_governc $cov i.year i.year#city_new if retire==0 , absorb(communityID) vce(cluster communityID)
psacalc delta anti_governc, rmax(0.341)
psacalc delta anti_governc, rmax(0.4)
psacalc delta anti_governc, rmax(0.5)

** Quantity of liquor
areg logliquor governc anti_governc $cov i.year i.year#city_new if retire==0 , absorb(communityID) vce(cluster communityID)
psacalc delta anti_governc, rmax(0.393)
psacalc delta anti_governc, rmax(0.4)
psacalc delta anti_governc, rmax(0.5)

//Table 8 //
use "$data_dir\data.dta",clear
set more off

global cov age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area

** Number of cifarettes
areg logcig1 governc anti_governc $cov i.year i.year#city_new if retire==0 , absorb(communityID) vce(cluster communityID)
psacalc beta anti_governc, rmax(0.409)
psacalc beta anti_governc, delta(0.9) rmax(0.409)
psacalc beta anti_governc, delta(0.8) rmax(0.409)
psacalc beta anti_governc, delta(0.7) rmax(0.409)
psacalc beta anti_governc, delta(0.6) rmax(0.409)
psacalc beta anti_governc, delta(0.5) rmax(0.409)
psacalc beta anti_governc, delta(0.4) rmax(0.409)
psacalc beta anti_governc, delta(0.3) rmax(0.409)
psacalc beta anti_governc, delta(0.2) rmax(0.409)
psacalc beta anti_governc, delta(0.1) rmax(0.409)
psacalc beta anti_governc, delta(0) rmax(0.409)

** Expenditure of cigarettes
areg logcig3 governc anti_governc $cov i.year i.year#city_new if retire==0 , absorb(communityID) vce(cluster communityID)
psacalc beta anti_governc, rmax(0.409)
psacalc beta anti_governc, delta(0.9) rmax(0.434)
psacalc beta anti_governc, delta(0.8) rmax(0.434)
psacalc beta anti_governc, delta(0.7) rmax(0.434)
psacalc beta anti_governc, delta(0.6) rmax(0.434)
psacalc beta anti_governc, delta(0.5) rmax(0.434)
psacalc beta anti_governc, delta(0.4) rmax(0.434)
psacalc beta anti_governc, delta(0.3) rmax(0.434)
psacalc beta anti_governc, delta(0.2) rmax(0.434)
psacalc beta anti_governc, delta(0.1) rmax(0.434)
psacalc beta anti_governc, delta(0) rmax(0.434)

** Frequency of liquor
areg f_liquor governc anti_governc $cov i.year i.year#city_new if retire==0 , absorb(communityID) vce(cluster communityID)
psacalc beta anti_governc, rmax(0.341)
psacalc beta anti_governc, delta(0.9) rmax(0.341)
psacalc beta anti_governc, delta(0.8) rmax(0.341)
psacalc beta anti_governc, delta(0.7) rmax(0.341)
psacalc beta anti_governc, delta(0.6) rmax(0.341)
psacalc beta anti_governc, delta(0.5) rmax(0.341)
psacalc beta anti_governc, delta(0.4) rmax(0.341)
psacalc beta anti_governc, delta(0.3) rmax(0.341)
psacalc beta anti_governc, delta(0.2) rmax(0.341)
psacalc beta anti_governc, delta(0.1) rmax(0.341)
psacalc beta anti_governc, delta(0) rmax(0.341)

** Quantity of liquor
areg logliquor governc anti_governc $cov i.year i.year#city_new if retire==0 , absorb(communityID) vce(cluster communityID)
psacalc beta anti_governc, rmax(0.393)
psacalc beta anti_governc, delta(0.9) rmax(0.393)
psacalc beta anti_governc, delta(0.8) rmax(0.393)
psacalc beta anti_governc, delta(0.7) rmax(0.393)
psacalc beta anti_governc, delta(0.6) rmax(0.393)
psacalc beta anti_governc, delta(0.5) rmax(0.393)
psacalc beta anti_governc, delta(0.4) rmax(0.393)
psacalc beta anti_governc, delta(0.3) rmax(0.393)
psacalc beta anti_governc, delta(0.2) rmax(0.393)
psacalc beta anti_governc, delta(0.1) rmax(0.393)
psacalc beta anti_governc, delta(0) rmax(0.393)

//Table 9 //
use "$data_dir\data.dta",clear
set more off

global cov age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area

** Panel A
	reghdfe logcig1 governc anti_governc $cov if retire==1,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m1
   reghdfe logcig3 governc anti_governc $cov if retire==1,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m2
   reghdfe f_liquor governc anti_governc $cov if retire==1,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m3
   reghdfe logliquor governc anti_governc $cov if retire==1,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m4
   
** Panel B
	reghdfe logcig1 governc official_governc $cov if retire==1,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m5
   reghdfe logcig3 governc official_governc $cov if retire==1,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m6
   reghdfe f_liquor governc official_governc $cov if retire==1,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m7
   reghdfe logliquor governc official_governc $cov if retire==1,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m8

** Panel C
	reghdfe logcig1 governc official_s_governc $cov if retire==1,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m9
   reghdfe logcig3 governc official_s_governc $cov if retire==1,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m10
   reghdfe f_liquor governc official_s_governc $cov if retire==1,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m11
   reghdfe logliquor governc official_s_governc $cov if retire==1,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m12
   
    esttab m1 m2 m3 m4 m5 m6 m7 m8 m9 m10 m11 m12 using "Table_9.csv",append order(governc anti_governc official_governc official_s_governc) legend starlevels( * 0.10 ** 0.05 *** 0.010) stats(N r2, star(r2) fmt(0 3) labels("Observations" "R-sq")) b(%9.3f) se(%9.3f) nogaps keep(governc anti_governc official_governc official_s_governc age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area _cons)	
	
//Table 10 //
use "$data_dir\data.dta",clear
set more off

global cov age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area
	
   reghdfe logcig1 firm anti_firm $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m1
   reghdfe logcig3 firm anti_firm $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m2
   reghdfe f_liquor firm anti_firm $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m3
   reghdfe logliquor firm anti_firm $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m4
    
   reghdfe logcig1 firm official_firm $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m5
   reghdfe logcig3 firm official_firm $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m6
   reghdfe f_liquor firm official_firm $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m7
   reghdfe logliquor firm official_firm $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m8  
   
   reghdfe logcig1 firm official_s_firm $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m9
   reghdfe logcig3 firm official_s_firm $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m10
   reghdfe f_liquor firm official_s_firm $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m11
   reghdfe logliquor firm official_s_firm $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
   est store m12
	
	esttab m1 m2 m3 m4 m5 m6 m7 m8 m9 m10 m11 m12 using "Table_10.csv",append order(firm anti_firm  official_firm official_s_firm) legend starlevels( * 0.10 ** 0.05 *** 0.010) stats(N r2, star(r2) fmt(0 3) labels("Observations" "R-sq")) b(%9.3f) se(%9.3f) nogaps keep(firm anti_firm  official_firm official_s_firm age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area _cons)
	
//Table 11 //
use "$data_dir\data.dta",clear
set more off

global cov age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area

** Panel A	
	reghdfe logcig1 governc anti_governc $cov if retire==0 & before_retire5_max == 1,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m1
	reghdfe logcig3 governc anti_governc $cov if retire==0 & before_retire5_max == 1,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m2
	reghdfe f_liquor governc anti_governc $cov if retire==0 & before_retire5_max == 1,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m3
	reghdfe logliquor governc anti_governc $cov if retire==0 & before_retire5_max == 1,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m4
	reghdfe logcig1 governc anti_governc $cov if retire==0 & before_retire5_max!= 1 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m5
	reghdfe logcig3 governc anti_governc $cov if retire==0 & before_retire5_max!= 1 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m6
	reghdfe f_liquor governc anti_governc $cov if retire==0 & before_retire5_max!= 1 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m7
	reghdfe logliquor governc anti_governc $cov if retire==0 & before_retire5_max!= 1 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m8
   
   esttab m1 m2 m3 m4 m5 m6 m7 m8 using "Table_11A.csv",append order(governc anti_governc) legend starlevels( * 0.10 ** 0.05 *** 0.010) stats(N r2, star(r2) fmt(0 3) labels("Observations" "R-sq")) b(%9.3f) se(%9.3f) nogaps keep(governc anti_governc age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area _cons)

** Panel B
   reghdfe logcig1 governc official_governc $cov if retire==0 & before_retire5_max == 1,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m1
	reghdfe logcig3 governc official_governc $cov if retire==0 & before_retire5_max == 1,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m2
	reghdfe f_liquor governc official_governc $cov if retire==0 & before_retire5_max == 1,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m3
	reghdfe logliquor governc official_governc $cov if retire==0 & before_retire5_max == 1,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m4
	reghdfe logcig1 governc official_governc $cov if retire==0 & before_retire5_max!= 1 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m5
	reghdfe logcig3 governc official_governc $cov if retire==0 & before_retire5_max!= 1 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m6
	reghdfe f_liquor governc official_governc $cov if retire==0 & before_retire5_max!= 1 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m7
	reghdfe logliquor governc official_governc $cov if retire==0 & before_retire5_max!= 1 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m8
   
   esttab m1 m2 m3 m4 m5 m6 m7 m8 using "Table_11B.csv",append order(governc official_governc) legend starlevels( * 0.10 ** 0.05 *** 0.010) stats(N r2, star(r2) fmt(0 3) labels("Observations" "R-sq")) b(%9.3f) se(%9.3f) nogaps keep(governc official_governc age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area _cons)

** Panel C
	reghdfe logcig1 governc official_s_governc $cov if retire==0 & before_retire5_max == 1,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m1
	reghdfe logcig3 governc official_s_governc $cov if retire==0 & before_retire5_max == 1,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m2
	reghdfe f_liquor governc official_s_governc $cov if retire==0 & before_retire5_max == 1,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m3
	reghdfe logliquor governc official_s_governc $cov if retire==0 & before_retire5_max == 1,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m4
	reghdfe logcig1 governc official_s_governc $cov if retire==0 & before_retire5_max!= 1 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m5
	reghdfe logcig3 governc official_s_governc $cov if retire==0 & before_retire5_max!= 1 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m6
	reghdfe f_liquor governc official_s_governc $cov if retire==0 & before_retire5_max!= 1 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m7
	reghdfe logliquor governc official_s_governc $cov if retire==0 & before_retire5_max!= 1 ,absorb(year communityID year#city_new) vce(cluster communityID)
	est store m8
   
   esttab m1 m2 m3 m4 m5 m6 m7 m8 using "Table_11C.csv",append order(governc official_s_governc) legend starlevels( * 0.10 ** 0.05 *** 0.010) stats(N r2, star(r2) fmt(0 3) labels("Observations" "R-sq")) b(%9.3f) se(%9.3f) nogaps keep(governc official_s_governc age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area _cons)

//Table 12 //
use "$data_dir\data.dta",clear
set more off
set seed 12345
set sortseed 12345
duplicates drop householdID year,force

tobit log_gift governc anti_governc $cov i.year i.communityID_new i.year#i.city_new if retire==0, ll(0) vce(cluster communityID)
margins, dydx(anti_governc)
margins, dydx(governc)

tobit log_gift governc official_governc $cov i.year i.communityID_new i.year#i.city_new if retire==0, ll(0) vce(cluster communityID)
margins, dydx(official_governc)
margins, dydx(governc)

tobit log_gift governc official_s_governc $cov i.year i.communityID_new i.year#i.city_new if retire==0, ll(0) vce(cluster communityID)
margins, dydx(official_s_governc)
margins, dydx(governc)
	
//Table 13 //
use "$data_dir\data.dta",clear
set more off

global cov age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area	
	
reghdfe friends governc anti_governc $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
est store m1
reghdfe majio governc anti_governc $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
est store m2
reghdfe shequ governc anti_governc $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
est store m3
reghdfe friends governc official_governc $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
est store m4
reghdfe majio governc official_governc $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
est store m5
reghdfe shequ governc official_governc $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
est store m6
reghdfe friends governc official_s_governc $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
est store m7
reghdfe majio governc official_s_governc $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
est store m8
reghdfe shequ governc official_s_governc $cov if retire==0,absorb(year communityID year#city_new) vce(cluster communityID)
est store m9

esttab m1 m2 m3 m4 m5 m6 m7 m8 m9 using "Table_13.csv",append order(governc anti_governc) legend starlevels( * 0.10 ** 0.05 *** 0.010) stats(N r2, star(r2) fmt(0 3) labels("Observations" "R-sq")) b(%9.3f) se(%9.3f) nogaps keep(governc anti_governc official_governc official_s_governc age age2 edu1 marrital gender log_wage_transfer fami_size area _cons)
	
//Figure A2 //	
use "$data_dir\data_PSM.dta",clear
set more off		

** The first graph	
permute antinew_governc beta = _b[antinew_governc] se = _se[antinew_governc] df = e(df_r), ///
 reps(500) rseed(123) saving("simulations_lnNcig.dta"): ///
  reghdfe logcig1 governc antinew_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)

use "simulations_lnNcig.dta", clear
gen t_value = beta/se
gen p_value = 2*ttail(df,abs(beta/se))
twoway(scatter p_value beta)(kdensity beta ,yaxis(2))

twoway (scatter p_value beta, msymbol(smcircle_hollow) mcolor(black))(kdensity beta ,yaxis(2) lp(solid)) , ///
       title("") ///
       xlabel(-0.2(0.1)0.2 -0.1 "-0.1" -0.2 "-0.2" 0.1 "0.1" 0.2 "0.2",format(%7.1f) angle(0)) ///
       ylabel(0(0.2)1, format(%7.1f) angle(0) nogrid axis(1)) ///
       ylabel(0(2)6, format(%7.1f) angle(0)  nogrid axis(2)) ///
       xtitle("Regression coefficient of LnNcig") ///
       ytitle("P value" ,orientation(vertical) axis(1)) ///
       ytitle("Kernel density" ,orientation(vertical)  axis(2)) ///
       xline(0, lwidth(0.2) lp(dash))  ///
       xline(-0.20, lwidth(0.3) lp(solid)) ///
       yline(0.1,lwidth(0.2) lp(dash)) ///
       legend(label(1 "P value") label( 2 "Kernel density")) ///
       plotregion(style(none)) /// 
       graphregion(color(white))	
	
** The second graph	
permute antinew_governc beta = _b[antinew_governc] se = _se[antinew_governc] df = e(df_r), ///
 reps(500) rseed(123) saving("simulations_lnEcig.dta"): ///
  reghdfe logcig3 governc antinew_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)

use "simulations_lnEcig.dta", clear
gen t_value = beta/se
gen p_value = 2*ttail(df,abs(beta/se))
twoway(scatter p_value beta)(kdensity beta ,yaxis(2))

twoway (scatter p_value beta, msymbol(smcircle_hollow) mcolor(black))(kdensity beta ,yaxis(2) lp(solid)) , ///
       title("") ///
       xlabel(-0.4(0.2)0.4 -0.4 "-0.4" -0.2 "-0.2" 0.2 "0.2" 0.4 "0.4",format(%7.1f) angle(0)) ///
       ylabel(0(0.2)1, format(%7.1f) angle(0) nogrid axis(1)) ///
       ylabel(0(1)4, format(%7.1f) angle(0)  nogrid axis(2)) ///
       xtitle("Regression coefficient of LnEcig") ///
       ytitle("P value" ,orientation(vertical) axis(1)) ///
       ytitle("Kernel density" ,orientation(vertical)  axis(2)) ///
       xline(0, lwidth(0.2) lp(dash))  ///
       xline(-0.324, lwidth(0.3) lp(solid)) ///
       yline(0.1,lwidth(0.2) lp(dash)) ///
       legend(label(1 "P value") label( 2 "Kernel density")) ///
       plotregion(style(none)) /// 
       graphregion(color(white))  

** The third graph
permute antinew_governc beta = _b[antinew_governc] se = _se[antinew_governc] df = e(df_r), ///
 reps(500) rseed(123) saving("simulations_f_liquor.dta"): ///
  reghdfe f_liquor governc antinew_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)

use "simulations_f_liquor.dta", clear
gen t_value = beta/se
gen p_value = 2*ttail(df,abs(beta/se))
twoway(scatter p_value beta)(kdensity beta ,yaxis(2))

twoway (scatter p_value beta, msymbol(smcircle_hollow) mcolor(black))(kdensity beta ,yaxis(2) lp(solid)) , ///
       title("") ///
       xlabel(-4(2)4 -4 "-4" -2 "-2" 2 "2" 4 "4",format(%7.1f) angle(0)) ///
       ylabel(0(0.2)1, format(%7.1f) angle(0) nogrid axis(1)) ///
       ylabel(0(0.1)0.4, format(%7.1f) angle(0)  nogrid axis(2)) ///
       xtitle("Regression coefficient of Fliquor") ///
       ytitle("P value" ,orientation(vertical) axis(1)) ///
       ytitle("Kernel density" ,orientation(vertical)  axis(2)) ///
       xline(0, lwidth(0.2) lp(dash))  ///
       xline(-2.069, lwidth(0.3) lp(solid)) ///
       yline(0.1,lwidth(0.2) lp(dash)) ///
       legend(label(1 "P value") label( 2 "Kernel density")) ///
       plotregion(style(none)) /// 
       graphregion(color(white))  


 
** The forth graph
permute antinew_governc beta = _b[antinew_governc] se = _se[antinew_governc] df = e(df_r), ///
 reps(500) rseed(123) saving("simulations_logliquor.dta"): ///
  reghdfe logliquor governc antinew_governc $cov if retire==0 [pw=w_ate],absorb(year communityID year#city_new) vce(cluster communityID)

use "simulations_logliquor.dta", clear
gen t_value = beta/se
gen p_value = 2*ttail(df,abs(beta/se))
twoway(scatter p_value beta)(kdensity beta ,yaxis(2))

twoway (scatter p_value beta, msymbol(smcircle_hollow) mcolor(black))(kdensity beta ,yaxis(2) lp(solid)) , ///
       title("") ///
       xlabel(-1(0.5)1 -1 "-1" -0.5 "-0.5" 0.5 "0.5" 1 "1",format(%7.1f) angle(0)) ///
       ylabel(0(0.2)1, format(%7.1f) angle(0) nogrid axis(1)) ///
       ylabel(0(0.5)2, format(%7.1f) angle(0)  nogrid axis(2)) ///
       xtitle("Regression coefficient of LnQliquor") ///
       ytitle("P value" ,orientation(vertical) axis(1)) ///
       ytitle("Kernel density" ,orientation(vertical)  axis(2)) ///
       xline(0, lwidth(0.2) lp(dash))  ///
       xline(-0.237, lwidth(0.3) lp(solid)) ///
       yline(0.1,lwidth(0.2) lp(dash)) ///
       legend(label(1 "P value") label( 2 "Kernel density")) ///
       plotregion(style(none)) /// 
       graphregion(color(white)) 
	
//Figure A3 //
use "$data_dir\data.dta",clear
set more off

** Figure 2A
label define govern_c 0 "Non-Official" 1 "Official" 
label value govern_c govern_c
graph bar cig1 if retire==1 , over(govern_c) over (year) asyvars  ytitle(`"{stSerif:number/day}"',height(6)) title(`"{stSerif:A: Average number of cigarettes consumed}"') ylabel(0(2)10)

** Figure 2B
gen year1 = year if year == 2011 | year == 2013 | year == 2015
graph bar cig3 if retire==1, over(govern_c) over (year1) asyvars  ytitle(`"{stSerif:yuan/month}"',height(6)) title(`"{stSerif:B: Expenditure on cigarettes}"') ylabel(0(20)100)

** Figure 2C
graph bar f_liquor if retire==1, over(govern_c) over (year) asyvars  ytitle(`"{stSerif:number/month}"',height(6)) title(`"{stSerif:C: Average frequency of liquor drank}"') ylabel(0(5)15)

** Figure 2D
graph bar liquor if retire==1, over(govern_c) over (year) asyvars  ytitle(`"{stSerif:ml/month}"',height(6)) title(`"{stSerif:D: Quantity of liquor drank}"') ylabel(0(500)1500)	
	
//Table A2 //
* Kernel matching for the Ncig

use "$data_dir\data.dta",clear
set more off
	
keep if year == 2011

set seed 10101
gen ranorder = runiform()
sort ranorder

psmatch2 govern_c edu1 marrital gender health_condition log_wage_transfer fami_size area i.communityID_new i.year#city_new, outcome(cig1) kernel ate ties logit common

pstest edu1 marrital gender health_condition log_wage_transfer fami_size area, both graph
	
keep if _support == 1
keep ID _pscore
encode ID,gen (ID_new)

duplicates drop ID ,force

save "$data_dir\id_Ncig.dta",replace

use "$data_dir\id_Ncig.dta",clear

merge 1:m ID using "$data_dir\data.dta"
keep if _merge == 3
drop _merge

save "$data_dir\PSM_Ncig.dta",replace
	
* Kernel matching for the Ecig
use "$data_dir\data.dta",clear
set more off
	
keep if year == 2011

set seed 10101
gen ranorder = runiform()
sort ranorder

psmatch2 govern_c edu1 marrital gender health_condition log_wage_transfer fami_size area i.communityID_new i.year#city_new, outcome(cig3) kernel ate ties logit common

pstest edu1 marrital gender health_condition log_wage_transfer fami_size area, both graph	
	
keep if _support == 1
keep ID _pscore
encode ID,gen (ID_new)

duplicates drop ID ,force

save "$data_dir\id_Ecig.dta",replace

use "$data_dir\id_Ecig.dta",clear

merge 1:m ID using "$data_dir\data.dta"
keep if _merge == 3
drop _merge

save "$data_dir\PSM_Ecig.dta",replace
	
* Kernel matching for the Fliquor
use "$data_dir\data.dta",clear
set more off
	
keep if year == 2011

set seed 10101
gen ranorder = runiform()
sort ranorder

psmatch2 govern_c edu1 marrital gender health_condition log_wage_transfer fami_size area i.communityID_new i.year#city_new, outcome(f_liquor) kernel ate ties logit common

pstest edu1 marrital gender health_condition log_wage_transfer fami_size area, both graph		

keep if _support == 1
keep ID _pscore
encode ID,gen (ID_new)

duplicates drop ID ,force

save "$data_dir\id_Fliqour.dta",replace

use "$data_dir\id_Fliqour.dta",clear

merge 1:m ID using "$data_dir\data.dta"
keep if _merge == 3
drop _merge

save "$data_dir\PSM_Fliqour.dta",replace
	
* Kernel matching for the Qliquor
use "$data_dir\data.dta",clear
set more off
	
keep if year == 2011

set seed 10101
gen ranorder = runiform()
sort ranorder

psmatch2 govern_c edu1 marrital gender health_condition log_wage_transfer fami_size area i.communityID_new i.year#city_new, outcome(liquor) kernel ate ties logit common

pstest edu1 marrital gender health_condition log_wage_transfer fami_size area, both graph		

keep if _support == 1
keep ID _pscore
encode ID,gen (ID_new)

duplicates drop ID ,force

save "$data_dir\id_Qliqour.dta",replace

use "$data_dir\id_Qliqour.dta",clear

merge 1:m ID using "$data_dir\data.dta"
keep if _merge == 3
drop _merge

save "$data_dir\PSM_Qliqour.dta",replace


//Table A3 //

** Ncig
use "$data_dir\PSM_Ncig.dta",clear

** For the Difference between treated and untreated sample
ttable2 cig1 cig3 f_liquor liquor anti official official_s governc age edu1 marrital gender health_condition wage_transfer fami_size area if retire == 0,by(governc)

ttest cig1 if retire == 0,by(governc)
ttest cig3 if retire == 0,by(governc)
ttest f_liquor if retire == 0,by(governc)
ttest liquor if retire == 0,by(governc)
ttest anti if retire == 0,by(governc)
ttest official if retire == 0,by(governc)
ttest official_s if retire == 0,by(governc)
ttest governc if retire == 0,by(governc)
ttest age if retire == 0,by(governc)
ttest edu1 if retire == 0,by(governc)
ttest marrital if retire == 0,by(governc)
ttest gender if retire == 0,by(governc)
ttest health_condition if retire == 0,by(governc)
ttest wage_transfer if retire == 0,by(governc)
ttest fami_size if retire == 0,by(governc)
ttest area if retire == 0,by(governc)

** Ecig
use "$data_dir\PSM_Ecig.dta",clear

** For the Difference between treated and untreated sample
ttable2 cig1 cig3 f_liquor liquor anti official official_s governc age edu1 marrital gender health_condition wage_transfer fami_size area if retire == 0,by(governc)

ttest cig1 if retire == 0,by(governc)
ttest cig3 if retire == 0,by(governc)
ttest f_liquor if retire == 0,by(governc)
ttest liquor if retire == 0,by(governc)
ttest anti if retire == 0,by(governc)
ttest official if retire == 0,by(governc)
ttest official_s if retire == 0,by(governc)
ttest governc if retire == 0,by(governc)
ttest age if retire == 0,by(governc)
ttest edu1 if retire == 0,by(governc)
ttest marrital if retire == 0,by(governc)
ttest gender if retire == 0,by(governc)
ttest health_condition if retire == 0,by(governc)
ttest wage_transfer if retire == 0,by(governc)
ttest fami_size if retire == 0,by(governc)
ttest area if retire == 0,by(governc)

//Table A4 //

** Fliqour
use "$data_dir\PSM_Fliqour.dta",clear

** For the Difference between treated and untreated sample
ttable2 cig1 cig3 f_liquor liquor anti official official_s governc age edu1 marrital gender health_condition wage_transfer fami_size area if retire == 0,by(governc)

ttest cig1 if retire == 0,by(governc)
ttest cig3 if retire == 0,by(governc)
ttest f_liquor if retire == 0,by(governc)
ttest liquor if retire == 0,by(governc)
ttest anti if retire == 0,by(governc)
ttest official if retire == 0,by(governc)
ttest official_s if retire == 0,by(governc)
ttest governc if retire == 0,by(governc)
ttest age if retire == 0,by(governc)
ttest edu1 if retire == 0,by(governc)
ttest marrital if retire == 0,by(governc)
ttest gender if retire == 0,by(governc)
ttest health_condition if retire == 0,by(governc)
ttest wage_transfer if retire == 0,by(governc)
ttest fami_size if retire == 0,by(governc)
ttest area if retire == 0,by(governc)

** Qliqour
use "$data_dir\PSM_Qliqour.dta",clear

** For the Difference between treated and untreated sample
ttable2 cig1 cig3 f_liquor liquor anti official official_s governc age edu1 marrital gender health_condition wage_transfer fami_size area if retire == 0,by(governc)

ttest cig1 if retire == 0,by(governc)
ttest cig3 if retire == 0,by(governc)
ttest f_liquor if retire == 0,by(governc)
ttest liquor if retire == 0,by(governc)
ttest anti if retire == 0,by(governc)
ttest official if retire == 0,by(governc)
ttest official_s if retire == 0,by(governc)
ttest governc if retire == 0,by(governc)
ttest age if retire == 0,by(governc)
ttest edu1 if retire == 0,by(governc)
ttest marrital if retire == 0,by(governc)
ttest gender if retire == 0,by(governc)
ttest health_condition if retire == 0,by(governc)
ttest wage_transfer if retire == 0,by(governc)
ttest fami_size if retire == 0,by(governc)
ttest area if retire == 0,by(governc)

//Table A5 //	
use "$data_dir\data_PSM.dta",clear
set more off

global cov age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area 

	reghdfe logcig1 governc anti_governc $cov if retire==0 [pw=w_ate], absorb(year communityID c.year#city_new) vce(cluster communityID)
	est store m1
	reghdfe logcig3 governc anti_governc $cov if retire==0 [pw=w_ate], absorb(year communityID c.year#city_new) vce(cluster communityID)
	est store m2
	reghdfe f_liquor governc anti_governc $cov if retire==0 [pw=w_ate], absorb(year communityID c.year#city_new) vce(cluster communityID)
	est store m3
	reghdfe logliquor governc anti_governc $cov if retire==0 [pw=w_ate], absorb(year communityID c.year#city_new) vce(cluster communityID)
	est store m4
	
	esttab m1 m2 m3 m4 using "Table_A5.csv",append order(governc anti_governc antinew_governc) legend starlevels( * 0.10 ** 0.05 *** 0.010) stats(N r2, star(r2) fmt(0 3) labels("Observations" "R-sq")) b(%9.3f) se(%9.3f) nogaps keep(governc anti_governc antinew_governc age age2 edu1 marrital gender health_condition log_wage_transfer fami_size area _cons)	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	


